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1. Introduction 

Particularly after the discovery of a Higgs-candidate resonance at ATLAS and CMS, the 
precise description of multi-jet Standard Model (SM) processes at the LHC remains crucial. 
Major progress has recently been made in combining calculations of next-to-leading order 
(NLO) perturbative QCD corrections with Parton Shower (PS) based Monte Carlo event 
generators [1—8]. 

In light of the rapid succession of publications on merging multiple NLO calculations 
with event generators [1 , 6-8] , it seems hardly arguable that this long-standing issue is re- 
solved. From the point of next-to-leading order accuracy, scepticism about state-of-the-art 
methods is baseless. It has however also been pointed out that sub-leading logarithmic 
enhancements could be left after the merging procedure [9-11]. It has to be stressed that 
initially, this problem is not caused by the extension of tree- level methods to NLO, but 
already appears for CKKW- inspired tree-level merging schemes [12-16]. The introduc- 
tion of a merging scale (t MS ) introduces logarithmic dependencies, L = In fj,p/t MS , with a 
dominant contribution a™L 2n , in the multi-jet tree- level configurations. These terms are 
partly cancelled by the parton shower higher-order corrections to low-multiplicity states. 
Let us look at one-jet merging in W-boson production, with a strictly leading- logarithmic 
parton shower. Integrating the W+jet matrix element (ME) over the one-jet phase space 
introduces the dependencies a s L 2 and a s L l , while unresolved PS emissions in the zero-jet 
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state produce a s L 2 terms, but with negative sign. In the total cross section, the a s L 2 
terms cancel, and only a t M s-dependence beyond the accuracy of the parton shower (i.e. 
a s L 1 -terms) remains. This example is rather academic, since modern parton showers in- 
clude a matrix-element correction for W+jet production, thus ensuring that in one-jet 
merging, all dependencies on t MS are cancelled 1 . After the inclusion of matrix element 
corrections [17-20], the merging scale dependence enters for two-jet merging. However, 
since control of beyond-leading logarithmic contributions cannot be universally exerted in 
a parton shower, it is clear that dependencies enter at some level, commonly hoped to 
be at next-to-next-to-leading logarithmic accuracy. The logarithmic structure of a two-jet 
tree-level cross section for example has only an approximate equivalent in parton showers, 
meaning that certainly at 0(a 2 ), sub- leading contributions are not fully cancelled. This is 
often, for lack of a better term, referred to as violation of PS unitarity. 

We believe it important to investigate this issue more closely. In this publication, we 
take a step back from the remarkable progress in NLO merging and outline a tree-level 
merging method that amends the behaviour of parton showers to ensure that no spurious 
logarithmic enhancements are introduced by including multi-jet matrix elements. The 
foundation of this method is PS unitarity, i.e. the requirement that the lowest order cross 
section remains unchanged by methods introduced to ameliorate the description of shapes 
of observable. It should be noted that in particular the GKS matrix-element-correction 
method in VlNCIA [21,22] has emphasised a unitarity-based approach before. The aim of 
this article is more modest, in that we offer a novel prescription of combining the input 
used in multi-jet merging procedures, moving from the additive scheme of CKKW-inspired 
methods to an add-and-subtract method that preserves the total inclusive cross section. 

This article is structured as follows. Section ||| is intended as introduction to our 
perspective on PS unitarity. This will be succeeded by a brief discussion of tree-level 
merging in the CKKW-L scheme [13, 14, 16] and its problems in section ||, before we move 
on to construct a LO merging method that preserves PS unitarity in section |j. Section 
|5| presents results for including additional jets in W-production and QCD dijet processes 
in the novel procedure, which we call UMEPS (Unitary Matrix Element + Parton Shower 
merging). Finally, we give a discussion in section ^ and conclude in section fj]. 

2. Parton shower unitarity 

Without any outside intervention, parton showers act on a lowest order seed cross section as 
a unitary operator. In other words, showering dresses the constituents of a perturbatively 
calculated 2—7-2 process with radiation, in order to set the stage for hadronisation. By 
generating soft and collinear emissions, parton showering sums (at least) leading double- 
logarithmic enhancements to all orders. 

Before detailing how this is achieved, let us introduce some notation in order to make 
the formulae less cluttered. We will also classify partons to be either resolved jets or 
unresolved jets at a particular scale p MS . By this, we mean that in a jet algorithm that 

1 Note that the introduction of ME corrections already endows the PS with the tree-level description of 
the hardest jet, so that one-jet merging would not be necessary. 
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exactly inverts the parton shower, a parton would be resolved as a jet if the evolution 
scale at which it was emitted (p), as well as jet separations that have been changed by the 
emission due to recoil effects, are above p MS . F° r the parton shower to be invertible in a 
well-defined way, we assume the existence of on-shell intermediate states between splittings. 
The index MS foreshadows the use of this jet definition as merging scale 2 . 

The parton shower approximates the effect of virtual corrections on observables sen- 
sitive to the scales pi and pi + \ by integrating DGLAP splitting kernels P (z) over the 
unresolved phase space. This gives a term 



"Pi 



CtGfoutgoing} J ttG{incoming} J i Z ^ 

dpdza s (p)V i+1 (z,p) , (2.1) 
p»+i 

where the first terms on the left-hand side sums all possible unresolved final state emis- 
sions, and the second term includes all unresolved initial state splittings. The ratios of 
parton distribution functions / is only absent if no initial parton is taking part in the 
(unresolved) emissions. The Vi+i (and P) notation is rather symbolic to permit a certain 
degree of simplicity. We include symmetry factors and the typical p- and z-fractions from 
approximating the matrix element or multiplying Jacobian factors in P and Vi+i, e.g. for 
an initial state splitting q — > qg, we would have P = ^ - ■ It can be shown, by perform- 
ing the z-integration for a specific splitting kernel P, that the PS "virtual corrections" in 



eq. (2.1) indeed capture the leading logarithmic contributions of virtual corrections. Let 



us introduce the short-hand 

fi(x i ,p i ) = f+(x+,p i )fr(x7,p i ) (2.2) 

The Sudakov form factor, resumming unresolved emissions between scales pi and pt+i, is 
given by 

fit^i Pi) 

A s+i (xi, pupi+i) = 1 , u — rlls (xi, Pi,p i+ i) (2.3) 

f i+ l(Xi, pi+l) 

where ILs +i (x{, pi,pi + \) is the probability of no emission from state S + i between the pi and 
Pi+i- The no-emission probability in turn can be expressed as 

n,s +1 (xi,pi,p i+ i) = exp I - j dpdza s (p) V i+1 (z, p) > (2.4) 

We have kept Xi as an argument to remember that n,s +i (xi, pi, pi+%) contains Xj-dependent 
PDF factors through V%-\-\- 



2 We here restrict ourselves to a particular jet (merging scale) definition, for the sake of clarity. We 
further assumed that the PS evolution variable is a measure of "hardness", i.e. that soft and collinear 
divergences are located at p — > 0. All the following arguments apply for a general merging scale tMS- 
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Let us consider the case when no PS emission above a scale p MS is generated. The 
parton shower approximation of the resummed exclusive zero-jet cross section is then 



dcr ex 

~ /oOo,Po) \M (p , p R )\ 2 U s+0 (x ,po,Pms) 



fo{xo,Po) \Mo (po,pr)\ 2 (2.5) 



,2 rP ° 



fo(xo,Po) \M (po,Pr)\ / dpdza s (p)T>i(z,p) 

J PMS 



PO 



+ fo(x Q ,po) \M (po,Hr)\ 2 I dp 1 dz 1 a s (pi)Vi(z 1 ,p 1 ) 



Pus 



rpi 

x / dp 2 dz2a s (p2)Vi{z2,P2) + C(«s) 

J PMS 



where we have used the fact that 



2 



1 / f Po \ 

I dpdza s (p)V 1 (z,p)\ (2.6) 

PMS / 

Po rpi 
dpidz 1 a s {p 1 )T>i{zi,pi) I dp 2 dz 2 a s (p 2 ) V\ {z 2 , P2) ■ 

PMS •'PMS 

The zero-jet PS cross section is exclusive in the sense that no resolved emissions (i.e. 
emissions above p MS ) are produced. Beside resumming unresolved contributions to the 
zero-jet cross section, the parton shower also produces resolved emissions. The parton 
shower approximation to the cross section for emitting the hardest jet at scale p\ is 

da in 

—j- = dp 1 dz 1 f (x ,po) \M (p ,PR)\ 2 a s (pi)Vi (zi,pi)Il s+0 (x ,Po,Pi) • 



V\ contains PDF ratios and kinematical factors (see definition 2.1). The cross section is 
labelled with in for "inclusive", because the emission of further jets below p%, but above 

da in 

p MS , is still allowed. Note that is also exclusive in the sense that no resolved emissions 
above pi - the scale of the first emission - are possible. In the following, we will always 
call a cross section inclusive if the parton shower can (at least in principle) produce further 
resolved emissions, and exclusive otherwise. 

If only zero- and one-jet states are generated, the total cross section is given by the 
sum of exclusive zero-jet and inclusive one-jet cross sections. Let us analyse the total cross 
section in the approximation of having maximally one parton shower emission. It useful to 
rewrite eq. (|2.5D with help of definition eq. (|2.4|): 



drr ex 

-f- = fo(x ,p )\M (po,VR)\ 2 (2-7) 

2 f p0 

- fo(xo,po) \Mo(po,Hr)\ / dpdza s (p)Vi(z,p)U s+0 {x,p ,p) 

J PMS 
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so that the total cross section is 



' rln ex f Hn™ 



J d<p (fo{x ,po)\Mo(po,VR)\' 



IH) 



fo(xo,Po) \M (po,p R )\ / dpdza s (p)Vi (z,p)U s (x,p Q ,p) 



+0 

PMS 
PO 

+0 



+ / d<f) fo(x , p Q ) \M (po,p R )\ 2 I dpdza s (p)T , (z,p)U s+0 (x ,po,p) 



PMS 

d<p fo{x ,po)\Mo(po,p R )\ 2 (2.8) 

Thus, if the parton shower would stop after the first emission, the total cross section is not 
changed by the application of the parton shower. 

This small sketch does not really ensure that the total cross section is preserved after 
PS resummation. Parton showering usually generates more than one emission, so that 
only being concerned with a single emission might not be enough. The above argument 
can however be extended to any number of emissions. As an example, assume the PS had 
generated two emissions. Then, the one-jet cross section becomes exclusive by demanding 
that only one resolved emission has been produced, 

—j- = dpidz 1 f (x ,p )\Mo(po,p R )\ 2 
d(f>o 

a s {fix)Vx (zi,pi)n 5+0 (x ,AbPi)n 5+1 (x 1 ,p 1 ,p MS ) (2.9) 

and we need to add the PS approximation to the two-jet cross section 
dcr in 

—j- = dpidz 1 dp 2 dz2fo(xo,po)\M (po,PR)\ 2 
dcpo 

x a s {pi)Vi {z u px)TL s+0 (x ,po,Px) 

x a s (p 2 )V 2 (z 2 ,p 2 )U s+1 (x 1 ,p 1 ,p 2 )@(px ~ Pi) (2.10) 



Now we rewrite |2.9| by expanding the second no-emission probability 

1 - dp 1 dz 1 fo(x ,p )\Mo(po,p R )\ 2 a s (p 1 )V 1 (z 1 ,p 1 )U s+0 (x ,po,pi) 



rpi 

1 - / dp 2 dz 2 a s (p 2 ) V 2 {z 2 ,p 2 ) 

J PMS 



rpi rpi 
+ / dp 2 dz 2 I dp 3 dz 3 a s (p 3 )'P 2 (z3,p 3 ) + 0(a s 

J OMR J PMS 



'PMS J PMS 

dpxdz 1 fo(x ,po) \M {po,p R )\ 2 a s {pi)Vi (z 1 ,p 1 )U s+0 (x ,po,Px) (2.11) 



rpi 

1- dpdza s (p)V 2 (z,p)U s+1 (xi,pi,p) 

J PMS 
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where we have again used definition 2A to derive the last equality. If parton showering 
stops after generating maximally two emissions, the total cross section is given by 



,»= i^(<£ + f*r. + f r*£) (2 . 12) 



By comparing the second term in eq. (EUD with eq. (ETlO|) , we see that any PS contribution 



of two resolved jets cancels with terms containing one resolved and one unresolved jet. 
Already earlier, we saw that contributions with one resolved jet cancel against terms with 
zero resolved jets. Thus, we again find 

#o/o(zo, Po) \M (p ,p R )\ 2 

It is easy to extend this argument to as many emissions as wanted: Whenever the parton 
shower generates one emission, the change in the total cross section is counteracted exactly 
by unresolved contributions to states with one emission less. There is no need to correct 
the PS approximation of the zero-jet exclusive cross section in the presence of two-jet states 
- the zero-jet resummation is oblivious of two-jet states. That parton showers are unitary 
is understandable directly from their construction, since the branching of an underlying 
n-jet state produces a n + 1-jet state, which overwrites (i.e. removes) the n-jet state in the 
wake of the branching. 



3. The problem with CKKW-L 

Let us now review tree-level matrix element merging, more specifically the CKKW-L 
scheme 3 . We here focus mainly on issues related to parton shower unitarity, and refer 
to [16] for a detailed description of CKKW-L in PYTHIA8 [23] and to [12-15] for a more 
general introduction. 

Matrix element merging procedures are designed to improve the PS description of 
multi-jet observables. For this purpose, tree-level matrix element (ME) calculations are 
combined with the parton shower, i.e. tree-level-weighted phase space points with m "hard 
process particles" and n additional partons are included in the shower. In the following, 
we will often use the terms state, event, configuration or the symbol S+ n to refer to these 
n + m— body phase space points. As can be inferred from the form of S+ n , we will most 
often understand that the state contains m hard process particles, but not mention these 
particles explicitly. 

A consistent merging removes all overlap between ME states and the PS approximation. 
This is ensured by introducing a phase space cut p MS to separate the ME region from the 
PS region, and applying no-emission probabilities. The cut dependence is minimised by 
weighting configurations above and below p MS in identical fashion. The CKKW-L scheme 
constructs and chooses a sequence of lower-multiplicity states (a so-called parton shower 
history) for each ME event, since factors need to be generated that would, in the parton 

3 Although most of what we discuss also applies to other CKKW-inspired merging schemes. 
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shower evolution, have contributed though intermediate stages. With the help of the 
history, ME events will be reweighted with 

xofo(x ,p ) ( tj %ifi( x i,P: 

W n — - ; r X 




\ Xi-lfi-l(Xi-l, Pi) J 

Yl^S +i -APi-l, Pi) ) X ^S +n {Pn, Pms) (3.1) 



8=1 



Xnfn(x n , pF 



n 



OtsiPi) Xi-lfi-l(Xi-l,Pi-l) . 

-U-S+i-AXi-i, Pi-i, Pi) 



a sME %i— lfi— l( x i— I; Pi 
x n s+n (x n , p n , Pms) , (3.2) 

where pi are the reconstructed splitting scales, and S + i the reconstructed intermediate 
states. The first PDF ratio in eq. (3.1) ensures that all ME configurations are normalised 
to the same total cross section, given by the lowest order Born-level matrix element. The 
PDF ratios in brackets account for PDF factors in the shower splitting probabilities Vi for 
initial state backward evolution. The running of a s is correctly included by the second 
bracket. Finally, double-counting is prevented by multiplying no-emission probabilities. 

Let us investigate how the CKKW-L merging prescription changes the lowest-order 
inclusive cross section. For simplicity, we will highlight merging matrix elements with up 
to two additional jets with parton showers. In the simplest conceivable case of one-jet 
merging, applying CKKW-L defines the cross sections 

da ex 

-jT- = fo(xo 1 po)\M (po,pii)\ 2 Ils + o{xo,Po,PMs) (3.3) 
= fo(xo,po) \Mo (po,Pr)\ 2 x (l- f dpdza s {p)Vi{z,p) + 0(a 2 s )\ 

\ "'PMS / 

dp 1 dz 1 f (x ,po) ^ svriy J1 )^ 1,riJ \Mi (po,Pr)\ z ^s +0 (x ,p ,pi) (3.4) 



— — = dpidzifo{xo,p )— -— r Ml PO.PB "b +0 

#0 a s {p R ) fo{x ,pi) 



It is crucial to note that the tree-level one-jet matrix element is in general different from the 
approximate PS splitting kernels. The inclusive lowest-order cross section is only preserved 
if 



2 f P ° 

fo(xo,Po) \Mo (po,Pr)\ / dpdza s (p)Vi{z,p) (3.5) 

-'Pms 

P °A A ft \ Mpl) fl( x ^Pl) | kA , N,2 
dpidzifo{x ,p ) -, r—, ? \Ml (PQ, fJLR)\ , 

' Pms a s {p R ) fo{xo,Pi) 

i.e. in the case where the first parton shower emission is distributed exactly according the 
one-jet matrix element. In this case, we would not have needed a merging prescription, 
since the PS would have already produced the correct result. 

Though correcting the first PS splitting to the full tree-level result is reasonably simple, 
correcting higher multiplicities requires significantly more work. The VlNCIA program aims 
at solving this issue [21,22]. In general however, we are currently forced to rely on tree-level 
merging to improve the descriptions of multi-jet observables. 
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If a first-splitting-corrected PS is available, unitarity violations will enter when includ- 
ing matrix elements for the next higher jet multiplicity. Since the case of two-jet merging 
is also instructive for later considerations, we will list the contributions to the cross section 
below. 



d<7 ex 

-jj- = fo(xo,Po)\Mo(po,HR)\ 2 U s+0 (xo,po,PMs) (3.6) 
= fo{xo,Po) \M (po,m)\ 2 (3-7) 

x ( 1 - / dpdza s (p) Vi(z,p) ILs+o (so, Po,p) 
\ J Pms 

da i X a a tf ,a s (p!) fi(xi,pi) 2 

— — = dpidzif (x , po)— -. — r-TT 7 \Mi {po,fJ,R)\ (3.8) 

dpo a s {PR) Jo{xo,pi) 

n s+0 (x ,p ,pi) U s+1 (x 1 ,p 1 ,p MS ) 

= dpidzif (x ,po) a ^ Pl \ ^^ Xl,Pl \ \Mi {p 0l p R )\ 2 Ils +0 (x ,p ,Pi) (3.9) 
a s {PR) Jo{xo,pi) 



x (l- / dpdza a (p)V 2 (z,p)U s+1 (xi.pi,^)) 

V J PMS / 



— — = dpidz 1 dp2dz 2 fo{x , po) — - — -— r — - — -— \M 2 {Po,Pr)\ (3.10) 

#o «s(w?) fo{xo,Pi) a s {p R ) fi(x 1 ,p 2 ) 

n 5+0 (x ,po,pi) U s+1 {xi,pi,p 2 ) 



For a first-splitting-corrected PS all contributions not containing II,sr +1 cancel between 3.7 



and 3.9, except for the lowest order inclusive cross section. Unitarity is then guaranteed if 



dpdz [ dp 2 dz 2 f (x ,po) as j Pl \ ^[ Xl,Pl \ \M\ (p , pr)\ 2 a s (p 2 )V 2 [z 2 ,p 2 ) 
pi J p MS a s {fj,R) fa{x ,pi) 

fis+o {xo, Po, p) n s+1 {xi,p, P2) 

P °a a f Pl A a tl \ a s(P) fi( x UP) a B (p 2 ) / 2 (^2,P2) , KA , \|2/ Q -m 

dpdz / dp 2 dz 2 f {xQ,p Q )— — -— -— — -— r \M 2 (po, m)\ ( 3 -H) 

Pl Jp MS a B (fj, R ) fo{x ,p) a s {fi R ) h{xi,P2) 

ris +0 (x , po, p) u s+1 (xi, p, p 2 ) 

For this, the splitting kernels need to exactly reproduce the matrix element, phase space 
must be fully covered by the parton shower, and the no-emission probabilities need to be 
produced identically in both cases. Particularly the requirement that the phase space is 
completely covered is problematic, since parton showers commonly fill only phase space 
regions in which consecutive emissions are ordered in a decreasing evolution variable. 

Clearly, eq. ( |3. 11 ) is not fulfilled in standard PS programs, which at best are correct 



to next-to-leading logarithmic (NLL) accuracy. This means that the dependence on the 
merging scale would vanish to order a^L 4 and a 2 L 3 , but that there will be a residual 
logarithmic dependence of order a 2 L 2 . 

In the next section, we would like to outline a method that sidesteps these problems 
by using multi-jet matrix elements from the very beginning to build the resummation for 
low-multiplicity states. 
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4. Concepts of UMEPS 



The main concept we would like to emphasise is that appropriately weighted matrix ele- 
ments with additional jets can be used to induce resummation in lower-multiplicity states. 

For example, one-jet inclusive cross sections (eq. fl2.7|) ) can, by integrating over the 
phase space of the resolved jet, be manipulated to induce resummation in zero-jet cross 
section. No parton shower resummation above p MS would then be necessary in zero-jet 
contributions. This means that we can reorder the parton shower formula for the inclusive 
cross section: 

Ida, 



d(poyfo(xo,po)\M (po,pB)\ 2 (4.1) 

2 f Po f da in " 

- fo{xo,po) \M {po,Pr)\ / dpdza s (p)V 1 (z,p)U s+0 (x ,po,p)+ 1 1 



PMS J a( P0 



and generate explicitly from by integrating over the emission phase space. When 

including one additional jet into the parton shower, we can explicitly preserve the inclusive 
cross section by adding the samples 



in 



7 

= fo(x ,Po)\M (p ,p R )\ 2 (4.2) 
d(p 

— — = dpidzi/o(x ,po) — / — s \Mi (po, hr)\ U s+0 {xo,po,Pi) (4.3) 

#o a s (p R ) Jo{xo,Pi) 

—rr l = -[ dpdzfo(x ,p ) ^P^ll^llfl \m x ( P0ifJlR )\ 2 n s (x ,p Q ,p) (4.4) 
#o J PMS a s (p R ) fo{x ,p) 

Before we continue, let us pause and investigate how we attach parton showers to these 
samples. In zero-jet contributions, the effect of parton showers above p MS is already in- 
cluded, so that we only need to start the parton shower at p MS . If the one-jet matrix element 
is the highest multiplicity sample, we allow the shower to generate emissions below pi, as 
in traditional merging. Since cross section changes from allowing e.g. two resolved jets 
cancel exactly with unresolved jets in the one-jet state (see eq. ( 2.10| ) and eq. ( [2.11 ) and 



the discussion following 2.12j ), allowing the shower to produce resolved emissions does not 
invalidate unitarity. 

We call this method UMEPS, for unitary matrix element + parton shower merging. 
In principle, this method is as easily generalisable as traditional merging techniques, and 
shows, on a more detailed level, difficulties reminiscent of CKKW-L. To particularise, let 
us have a look at how two additional jets can be included by UMEPS. Naively, we would 
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simply add 

dcrT j j j j ft ,a s {pi) h{xx,Pi) a s {p 2 ) f 2 {x 2 ,p 2 ) , s,2 ,„ - 

— f- = dp 1 dz 1 dp 2 dz2fo{x ,po) — 7 — r-777 r — 7 — VT7 V ^2 (po, Pi?) ( 4 -5) 

#0 «s(Pr) 7o(zo,Pi) <Wii) h{xi,P2) 

xIi s+0 (x ,po,pi)U s+1 {x 1 ,p 1 ,p 2 ) 

= -dpidzi \ dpdzf (x ,po) — 7 VTT \ — 7 TT7 v -^2 (po, Pi?) ( 4 - 6 ) 

x~tts +0 (zo,Po,Pi)n s+1 (xi,pi,p) 



and treat 4.5 as highest multiplicity sample. It is however possible that due to undoing 



recoil effects, states with jets below p MS are produced by performing the integration in 4J3 
In this case, we take these contributions to be corrections to the zero-jet cross section, and 
integrate twice. After this amendment, two-jet UMEPS merging contains the contributions 

dn in 

-f- = f (x ,po)\Mo(p ,p R )\ 2 (4.7) 

d ^j- = dpidzif (x , po) ° s< ] Pl \ j-p^r \Mi (po,pn)\ 2 ~ns +0 (x ,po,pi) (4.8) 
d<Po ® s {pr) j {x , pi) 

= - [ PO dpdzfo(xo,p ) a f P) J f l f l,P \ |-Mi(po,Pfi)| 2 n s+0 (xo,po,p) (4.9) 
d(po J PMS «s(Pi?) Mxq,p) 

da 2 n A a a a ti \ a s(Pi) hixupx) a s {p 2 ) f 2 (x 2 ,p 2 ) 2 (A ln x 

— — = dp 1 dzidp2dz2fo{xo,p )— — -— r— — r— r \M 2 [po, Pr)\ (4.10) 

dcpo <MPi?) fo(x ,px) a s (p R ) /i(x 1 ,p 2 ) 

Hs +0 (x ,P0,Pl)IIs +1 (Xi,pi,p2) 

<fr£»i , , [ P1 a a ft ^ Qs(pi) /i(a:i,pi) Qs(p) f2(x 2 ,p) 2 

— — — = -dp x dz x \ dpdzf (x ,p )— — -— - — — -— - |.M 2 (po, Pi?)| 

dcpo J PMS a s {p R ) fo(x ,pi) a s (p R ) h{xi,p) 

@(pi -p M s)n5 +0 (^o,po,pi)n s+1 (x!,pi,p) (4.11) 

da % 2 %Q f P0 , , f Pl » « s (Pa) /i(ari, p a ) a s (p b ) f 2 (x 2 ,p b ) 2 

j 2 = ~ / dpadz a / dp b dz b fo{x ,p )— r— r— r— r | M 2 {p , p R ) \ 

d(po J Pl J PMS a s (Pi?) 7o {%o, Pa) « S (P«) fi{xi, Pb) 

@(Pms - Pa)n 5+0 (x , p , p a ) U S+1 (Xi, p a ,Pb) (4.12) 

UMEPS can then be extended to arbitrary jet multiplicity. The main idea is that in order 
keep unitarity, we have to subtract all contributions that we add as higher multiplicity 
matrix elements. The subtractions are constructed with PS unitarity as a guideline. For 
brevity, we introduce the short-hands 

^f- = B n w' n = B n and = - / d n ' m ^ B n w' n = - / B„ 

dcpo dcpo J Js 

where w' n will be defined below. The symbol J s B„^ m indicates that more than one inte- 
grations had to be performed since all of the states S* +n _i, . . . , S +m+ i contained partons 
below the merging scale. The integration(s) will be achieved by substituting the input 
event with a reconstructed lower-multiplicity event of the parton shower history, as will be 



discussed in section 4.1. This substitution method is indicated by the subscript s on the 



integral sign. The weight w' n that needs to be applied to tree-level events to produce the 
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-TT- sample is given by 



/ %nfn{%ni Pn) 
W„ = — r X 



%nfri{%ni Pf) 



n 

i=i 



a s (pn) Xi-xfi-x(xi-\, pi) 



fts+i-i(xi-l, Pi-l, Pi) 



(4.13) 



j£n fn {%n ipn) %n fn (*£n iPfi) 

n 



iWi \ x i-\, Pi-\i Pi) 



i=l 

This weight differs from the CKKW-L weight in eq. ( |3.2| ), since it does not contain the 
last no-emission probability IIs +ti (i n , p n , Pms)> the last line in eq. ( |3 . 2| ) . In the UMEPS 
procedure, this factor is instead included by subtracting the integrated, reweighted, next- 
higher multiplicity sample, thus conserving unitarity in a way reminiscent of standard 
parton showers. The probability of having no resolved emissions off the zero-jet states in 
eq. ( |4.7[ ) for example, is included through the contributions in eqs. (PI) and eq. (|Og ). 

Armed with this notation, the prediction of an observable O in 2— jet merged UMEPS 
becomes 



(O) 



6o|o(5 + oi)B - O(S +0j ) Jb^o - O(S +0j ) ^ 
+ f 0(S +lj )B L - J 0(S +lj ) Jb^i 
+ II 0(S +2j )B 2 \ , 



B 



2^0 



(4.14) 



where we have used the notation S+ n j to indicate states with n resolved jets, resolved 
meaning above the cut p MS as defined by the merging scale definition. More generally, the 
outcome of merging n additional partons with the UMEPS method is 



(O) 



b yD(S +0 j 
+ 1 0(S +l3 ) 



Bo - / Bi-mi 

s 



B 



2^0 



Bi 



B 



2->l 



• • • — / Btv^o 
Bjv^i 



+ 



+ 



+ 



0(S +N - 



B/v- 



N-l 



B 



N->N-1 



0(S +Nj ) B 



N 

E 

n=0 



0(S, 



II J I 



N 



B n 



N 

E 

i=n+l 



B; 



(4.15) 



The generation of B n — and J s B n _j. m — events will be summarised in section LI. It should be 
noted that the treatment p MS — unordered integration results is heavily influenced by how 



- 11 - 



CKKW-L includes states with /9 MS — unordered emissions, which was discussed in detail 
in [16]. Precisely for states which evolve from a state below p MS to a state above p MS do 
CKKW-L and the truncated-shower [24] approach differ. It can thus be imagined that 
other ways of treating such notorious configurations show improved behaviour. For now, 
we will not discuss such possibilities, and instead, if necessary, integrate multiple times, 
until a state above p M s is produced. 




W" 



Figure 1: An example of a matrix element contribution without a complete shower history. If the 
parton shower does not include W— boson radiation, only the two gluon emissions can be reclustered, 
and cc — > udW~ has to be considered a separate hard process. 

A well-known challenge of merging prescriptions is the treatment of configurations that 
could never have been produced by a sequence of shower splittings. This can happen if 
the PS does not include all possible splittings of the model. Figure [| for example cannot 
be produced by a shower that does not allow W— boson radiation. Such states cannot 
be projected onto a lower-multiplicity underlying process, and will thus introduce small 
unitarity violations in UMEPS. 

4.1 Procedure step- by-step 

To implement UMEPS, we need to be able to perform the necessary integrations. Although 
the formulae could convey a feeling of complexity, these integrations are factually already 
needed in traditional merging approaches. All modern CKKW-inspired schemes need to 
construct a history of parton shower states for input matrix element events, because oth- 
erwise, no trial showers can be used to generate Sudakov form factors dynamically. Thus, 
a sequence of states S+ n — > S+ n -i —>••••—> 5+i — > S+o is always available. The sequence 
is constructed by inverting the shower mapping of radiative phase space on each state, i.e. 
a parton in S+ n is removed, and its momentum distributed amongst the remaining parti- 
cles, leading to a state S +n _i 4 . This is exactly the integration we need: To produce the 
integrated version of S +n , we simply replace it by 5 +n _i, but keep the full weight, k— fold 
integrations can be achieved by replacing S +n with 5 , +n _/ c . 

With this, we have all ingredients to construct the UMEPS scheme. If not mentioned 
explicitly, all weights in UMEPS are generated precisely as in CKKW-L. To avoid unnec- 

4 The mapping used for the current paper is given in Appendix B.2 of [16] 
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essary complications, we will here gloss over many technicalities that have already been 
addressed in CKKW-L, and are directly borrowed by UMEPS. A discussion of Sudakov 
reweighting for states without parton-shower ordered histories, for example, can be found 
in [16]. 

The UMEPS algorithm has two parts - a part in which we keep the matrix el- 
ement configurations (B n — events) and another in which we integrate over emissions 
(J s B n _> m -events): 

I. Produce Les Houches event files (LHEF) [25] with a matrix element generator for 
n = 0,1 ... N extra jets with a regularisation cut-off, p MS , typically using a fixed 
factorisation scale, pf, and a fixed a s (pR). 

II. Pick a jet multiplicity, n, and a state S +n according to the cross sections given by 
the matrix element generator. 

1. Find all shower histories for the state S +n , pick a sequence according to the 
product of splitting probabilities. Only pick un-ordered sequences if no ordered 
sequence was found. Only pick incomplete paths if no complete path was con- 
structed. 

2. Perform reweighting: For each ^ i < n, 

i. Start the shower off the state S + i at pi, generate a trial state R with scale 
Pr. If pn > Pi+i, veto the event and start again from II. 

ii. Calculate the weight factor 

_ a s (p i+ i) xff+{xf,pj) xlfr{xT,pi) 
a s {PR) x\ fi (x- ,p i+ i)x i f i (x- ,p i+l ) 

3. Start the shower from S+ n . 

i. If n < N, start the shower at p n , veto any shower emission producing an 
additional resolved jet. 

ii. If n = N, start the shower at p n . 

III. If the event was not rejected, multiply the event weight by 

/ _ x n fn i x n > Pn) x n fn i x n i Pn) TT , . -, 7 \ 

w n — +7+7 + \ X - f- ( - s x {^■ Ll ) 
Xn ]n \Xn ; PF ) X n J n \X n , Pf ) 

V. Start again from II. 

The second part, i.e. producing J g B n _^ m — events to effect lower-multiplicity PS resumma- 
tion, requires only two changes: 

II. 3 Replace the matrix-element state by S +n -i, or the first state S + i with all / < n — 1 
partons above the merging scale. If no integrated state can be constructed, i.e. if only 
incomplete paths were found, reject the event. For valid events, start the shower at 
p n , veto any shower emission producing an additional resolved jet. 
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III. If the event was not rejected, multiply the event weight by 

/ x nfn( x n'Pn) x nfn ( x n>Pn) TT /. -. q n 

~ W n = + ,+ / + T X \ X 11 W i ( 4 - 18 ) 

x n In \ x n ; PF ) x n In \ x n j PF ) 

Finally, all samples generated in the first and second parts are added to give the UMEPS 

dcr? x 

prediction. Note that in order produce correctly normalised cross sections we need 

to include ratios of parton distributions and a s ratios into the weight. This is analogous 



to the CKKW-L method (see the a s - and PDF-factors in eq. ( |3.2D ). It is worthwhile 
to notice that the UMEPS scheme can in principle be implemented by using structures 
already existing in traditional merging codes. Basically, compared to traditional merging, 
the B n — contributions do not carry a no-emission probability for emissions off the ME event. 
The J B„_> m - samples can easily be extracted from merging codes. 

How multiple partonic interactions (MPI) are added to the merged samples requires 
a short discussion. In principle, we stay true to the philosophy of the algorithm outlined 
in [16], i.e. we want to make sure that the merging method does not artificially suppress 
hard secondary scatterings, which in PYTHIA8 are interleaved with the parton shower. 
The interleaving means that the PS is competing with the MPI's, and the probability of 
emitting a parton in the PS is not only governed by the standard no-emission probability, 
n,s +n but is also multiplied by a no- MPI probability, Uf p *. Hence all m-jet (both B m and 
f B n ^. m ) samples generated by our algorithm above need to be multiplied by the no-MPI 
probabilities 

m— 1 

ncift.mi). (4.19) 

j=0 

which are easily incorporated in the trial showers described above. We also need to include 
the actual MPI's. Here the philosophy is that as soon as we have a MPI at some scale, we 
ignore corrections from the full tree-level matrix element on softer jets from the primary 
interaction, and allow them to be described by the PS alone. Hence, when we start the 
shower from a given m-parton state (with m < N) in step II. 3, we choose the reconstructed 
p m as starting scale. As described before, we veto any parton emission above p MS - However, 
if a MPI is generated above p MS , it is accepted and the shower is allowed to continue without 
any further veto. (For the m = N case, the shower including MPI is allowed without 
restrictions, starting from pn-) In this way we achieve the same goal as in [16]: If the 
n < N hardest jets in an event all belong to the primary interaction, they are described by 
the tree-level ME, while all other jets are given by the (interleaved) PS. Just as in [16], the 
treatment of pure QCD jet production means that the Born-level cross section is properly 
eikonalized by the no-MPI factor, by allowing MPI's all the way from y/s in the trial shower 

for n^pi. 

5. Results 

We have implemented UMEPS merging in PYTHIA8, and will make the necessary code 
public in the next major release version. In this section, we will concentrate on predictions 
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for W-boson and QCD jet production at the LHC. However, the code aims to achieve the 
same generality as the implementation of CKKW-L in PYTHIA8. 

All input matrix element configurations are taken from Les Houches Event Files gen- 
erated with MadGraph/MadEvent, with the following settings: 

• Fixed renormalisation scale ^ir = M|, fixed factorisation scale fj,p = for W- 
production. For 2—^2 processes in pure QCD, we use fi r ,2^2 = m± t im± > i and 
M/,2->2 = min{ml tl ,ml A }. 

• CTEQ6L1 parton distributions and a s (M 2 ) = 0.130. 

• The merging scale p MS is defined by the minimal PYTHIA8 evolution p±^jk of all 
possible combinations of three partons in the event. p±^jk for a single combination 
of three particles i, j and k is defined as 

Z ijk (l ~ Zijk)Qij With Qfj = {Pi+Pj) 2 , Zijk = x J k +x j<ik , 

Xljk = YTT'^ fo^FSR 



2 

P±,ijk 



II - z ijk )Ql with Q% = -(» - p 3 f , z tjk = { XX"r 
for ISR 



(5.1) 



• In QCD 2—7-2 scatterings, the kinematical transverse momentum of jets is required 
to be larger than pxj = 5 GeV. 

The value of a s (M 2 z ) was set to match the a s -value obtained in fitting the PDFs used in the 
ME calculation. To generate results, we have chosen the merging scale definition to closely 
match the parton shower evolution variable. The algorithm does however not depend on 
this particular choice. All jets needed for analysis purposes were defined with help of 
fast jet-routines [26]. The momentum of the intermediate W-boson will, if required, be 
extracted directly from the Monte Carlo event. We will compare UMEPS to the CKKW-L 
implementation in PYTHIA8. The problems we choose to highlight should be regarded as 
criticism of the implementation in PYTHIA8, rather than an assessment of CKKW-inspired 
methods in general. 

5.1 W-boson production 

We begin by comparing the result of the removal of a jet by integration with the corre- 
sponding parton shower contribution. This is useful to assess if performing the integration 
by the replacement <S+ n +i — >• S +n produces the desired results. 

In the left panel of Figure ^, we compare the integrated one-jet matrix element (i.e. the 
O [oq (/^))-term of eq. (|4.9|)) with the shower approximation of the O (a\ (/^))-term in 
zero-jet events. The second term is of course just the O (a* (/J-r)) -contribution in eq. ([O]). 
The rapidity of the W-boson is identical in these two samples because PYTHIA8 is already 
matrix-element corrected for W+ j-states. This demonstrates that in W-boson production, 
generating the no-emission probability in zero-jet states with PYTHIA8, or by a reweighted, 
integrated one-jet matrix element are both legitimate ways to produce the same factor. 
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Figure 2: Comparison between C(a s )-terms of the parton shower with integrated matrix ele- 
ments, for W-boson production in pp collisions at Eqm — 7000 GeV. The merging scale value 
is Pms = 15 GeV. Left panel: Rapidity of the W-boson, intended for comparison between the 
integrated one-jet matrix element (labelled [dai^a/dy w ] ol ^ al - ) ) and the O(o!g)-term of the no- 
emission probability for having no emission above pms radiated off zero-jet states above pms (labelled 
[dao/dy w ] ( a ij ). Right panel: Transverse momentum of the parton, for W + j production in pp 
collisions at Ecm — 7000 GeV, intended for comparison between the integrated two-jet matrix 
element (labelled [d<T2^>i / dp±] ( a2 - ) ) and the ©(a^-term of the no-emission probability for having 
no emission above p MS radiated off one-jet states, multiplying the W + j matrix element (labelled 
[dai/dpx] 0{aT) ). 



The right panel of Figure ^ investigates the difference between the parton shower 
approximation of no resolved emissions in one-jet states between the scales p\ and p M s and 
the result of constructing an unresolved emission by integrating over one parton in a two- 
jet matrix element. This means that we compare the one-jet matrix element, multiplied 
by the O (a\ (//R))-term of the no-emission probability Us +1 (xi, Pi, Pms) in eq. (3J3), with 



the O (ai? (//i?))-contribution in eq. ( pll ). The comparison shows that, as expected, the 
parton shower underestimates the hardness of the unresolved (second) emission, which is 
reminiscent of the fact that the inclusion of two-jet matrix elements into the PS prediction 
does in general increase the tail of the p± of the hardest jet. 

In Figure, |3] we show how matrix element samples contribute to this increase. All jet 
multiplicities enter, because the merging scale is not defined as the jet-separation of the 
/c_i_-algorithm, and since the merging scale cut acts on the matrix element state, while the 
jets are constructed from outgoing particles after the parton shower cascade. In CKKW-L, 
the high-p^ tail is dominated by the two-jet matrix element, with a major contribution 
from the one-jet states. The latter is significantly lower in UMEPS, a fact that we think 
crucial. UMEPS correctly cancels the inclusion of phase space points with two resolved jets 




Figure 3: Transverse momentum of the hardest jet, for W-boson production in pp collisions at 
Ecm = 7000 GcV, when merging up to two additional partons. Jets were defined with the kj_- 
algorithm, with k±^ m i n = 10 GcV. Multi-parton interactions and hadronisation were excluded. Left 
panel: Results of the CKKW-L scheme. The contributions are labelled Bo^o, Bi^i and B 2 W2 for 
CKKW-L-reweighted zero-, one- and two-jet matrix elements, respectively. Right panel: Results of 
the UMEPS scheme. The contributions are labelled Bo, Bi and B2 for UMEPS-reweighted zero-, 
one- and two-jet matrix elements, and / Bi^o and / B2->i for UMEPS-reweighted, integrated one- 
and two-jet samples. J s B2_j.o indicates the two-jet contribution that was integrated twice because 
the state 5+i after the first integration contained an unresolved jet. 

by using the two-jet matrix element to construct a better approximation of radiating an 
unresolved parton from one-jet states. We see both in CKKW-L (left panel of Figure ||), and 
in Figure ||, that the parton shower underestimates the hardness of two-parton states. The 
description of unresolved emissions enters into the no-emission probabilities, with a negative 
C(ag)-term. Thus, the contribution of showered one-jet states to the tail of p±i will be 
larger if the shower description of two-jet states underestimates hardness. UMEPS improves 
the description of the no-emission probability by ensuring that in inclusive observables, 
resolved two-jet states are cancelled, a feature that is at work in the tail of p±±. 

Variations in the description of p±x are also visible in Figure where we show the 
transverse momentum of the hardest jet in CKKW-L and UMEPS 5 . The trend sketched in 
the previous paragraph is particularly clear in the insets comparing to default PYTHIA8: 
UMEPS produces a softer tail in p L1 than CKKW-L. The harder tail in CKKW-L is 
due to a worse description of unresolved emissions. It is fair to say that the difference 
between CKKW-L and UMEPS hints at the size of relic effects from not cancelling the 
higher-multiplicity matrix elements in a well-defined way. Merging scale variations in tree- 
level merging schemes arise from a mismatch of unresolved emissions exponentiated in 
no-emission probabilities and tree-level matrix elements for hard, resolved jets. UMEPS 

Note that the co-variation of merged results in the ratio inset is due to fluctuations in the PYTHIA8 
reference curve. 
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Figure 4: Transverse momentum of the hardest jet, for W-boson production in pp collisions at 
Ecm = 7000 GeV, when merging up to two additional partons. Jets were defined with the fcj_- 
algorithm, with fcj_ !m i„ = 10 GeV. Multi-parton interactions and hadronisation were excluded. 
The lower insets show the deviation of merged results from default PYTHIA8, for three different 
Pivis-values. Left panel: Results of the CKKW-L scheme. Right panel: Results of the UMEPS 
scheme. 

has a significantly lower merging scale variation since the method enforces a cancellation 
of resolved and unresolved contributions. 

Figure |5| shows that for very exclusive observables, CKKW-L and UMEPS are virtually 
indistinguishable. In this example, this is of course expected since the treatment of the 
highest multiplicity (here, the two-jet) matrix element is identical for both cases. 

We would now like to perform a stress-test of the merging scale dependence. Since 
UMEPS properly cancels the effects of adding multi-jet matrix elements by subtracting 
their integrated counter-parts, it is in principle possible to push the merging scale to very 
small values. The variation of the inclusive cross section is shown in Figure ^. It is clear 
that UMEPS does indeed preserve the inclusive cross section, while for CKKW-L, very 
small merging scales lead to large changes, rendering the method unreliable. However, the 
error convergence in UMEPS is, due to the negative weights, significantly slower. We will 
comment on this below. 

These unitarity violations might not induce drastic effects in the description of hard- 
scale observables like the transverse-momentum distribution of the W-boson. However, 
magnifying the low-scale description of this observable (Figure |7|) reveals problems. Figure 
serves two purposes. It clearly shows that by pushing the merging scale to small values 
in CKKW-L, sub-leading contributions in the multi-jet matrix elements start to contribute 
more. Since those sub-leading contributions cannot be cancelled by the default parton 
shower, major increases over PYTHIA8 are found. UMEPS explicitly cancels these sub- 
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Figure 5: Transverse momentum of the second-hardest jet, for W-boson production in pp collisions 
at Ecm = 7000 GeV. Jets were defined with the fc_L-algorithm, with k± tm i n = 10 GeV, when merging 
up to two additional partons. Multi-parton interactions and hadronisation were excluded. The lower 
insets show the deviation of merged results from default PYTHIA8, for three different pMS-values. 
Left panel: Results of the CKKW-L scheme. Right panel: Results of the UMEPS scheme. 

leading terms, and thus leads to a reliable prediction. The second observation in Figure ^ 
is the dependence on the primordial transverse momentum parameter k± iP . This parameter 
was introduced in event generators to account for the transverse momentum of partons in 
the incoming protons, which cannot be generated in initial state DGLAP evolution. If this 
were the only effect to be modelled by k± tP , a value of k_\_. p « 0.3 GeV per incoming parton 
would seem appropriate. However, in current event generator tunes, significantly higher 
values « 2 GeV) are required [20,27], potentially to compensate for an incomplete 

phase space coverage in initial state showers due to the shower cut-off. The value of k± tP is 
mainly fixed by tuning to the position of the peak of the transverse momentum spectrum of 
the Z- or W boson. Increasing the value of kj_ tP roughly corresponds to pushing the peak 
to higher p± values. Figure |7| compares the UMEPS and CKKW-L predictions for the 
transverse momentum of the W boson, with k±„ p = 0.5 GeV, to default, tuned PYTHIA8 
with k± iP = 2.0 GeV. Unitarity violations in CKKW-L pull the peak back to lower p±. 
This fact is virtually unchanged if we had used /cj_ iP = 2.0 for CKKW-L predictions instead, 
suggesting that if we positively wanted to use a very low merging scale, an increase in 
kj_ tP would be necessary. UMEPS on the other hand can be used with very low merging 
scales, and in particular shows the interesting feature of matching the default PYTHIA8 
curve without having a high k± tP value. We believe this is due to a better modelling of 
logarithms of the form ln(l/x), which are present in the matrix element, and which are 
included in a unitary way in UMEPS - allowing for a much more natural value of k±^ p . 
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Figure 6: Inclusive cross section for UMEPS- and CKKW-L merging of up to two additional 
jets in W-boson production in pp collisions at Eqm = 7000 GeV (labelled <x mer gcd), m comparison 
to the the lowest-multiplicity inclusive cross section <7i nc i US i VC . The error bars represent only the 
statistical error on the merged cross section. For the UMEPS sample this becomes large for small 
merging scales, but as it is the same samples which are added and the subtracted, the central value 
stays very close to unity. 



This result is of course very preliminary, since there are e.g. correlations of the shower 
cut-off p± m m an d k±^ p . One would hope that matrix-element merging would allow to lower 
P±mim which might mean having to make a compromise for the value of k±^ p . We will 
come back to these aspects when presenting tunes for matrix-element-merged PYTHIA8 in 
a future publication. 



5.2 Dijet production 

We would further like to mention QCD dijet production at the LHC, since potential merging 
scale dependencies enter already when merging dijet- and three-jet matrix elements, and to 
demonstrate the flexibility of our implementation. The main objective of including QCD 
dijet production in this publication was to assess the treatment of MPI discussed at the 



end of section |4.1| . This is most effectively done by comparing to data, and before these 
we would like to only stress one issue. 

Figure |8] shows the transverse momentum of the third-hardest jet. We see that UMEPS 
and CKKW-L show very similar changes when compared to PYTHIA8. However, neither 
curves show the high-p^ increase seen in [16]. This is simply because we have revised the 
choice of the renormalisation scale in the core 2 — > 2 QCD scattering. In [16], the two 
powers of a s (fXu) in the core 2—^2 process were never touched, and thus were evaluated 
with the rather unfortunate choice hr = in the input Les Houches events. This 
has been rectified in the current release of (CKKW-L in) PYTHIA8, i.e. the scale choice 
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Figure 7: Transverse momentum of the W-boson, for W-boson production in pp collisions at 
Ecm = 7000 GeV, when merging up to two additional partons. Multi-parton interactions and 
hadronisation were excluded. The PYTHIA8 results are generated with the default settings, in 
particular with a primordial transverse momentum of k±. p = 2 GeV. All merged curves have been 
generated with k±. p = 0.5 GeV. The lower insets show the deviation of merged results from default 
PYTHIA8, for three different p MS -values. Left panel: Results of the CKKW-L scheme. Right panel: 
Results of the UMEPS scheme. 

(fi r ,2->2 = ™i,imi,i) is now dynamical. Potential reweighting (due to the usage of fixed 
in the LHEF generation is handled internally in PYTHIA8. The trend that pure QCD 
multi-jet matrix elements have a softer spectrum of well-separated jets has already been 
observed in [16,28]. The merging scale variation in UMEPS is within the statistical error 
of the samples. The statical uncertainty is larger in UMEPS than in CKKW-L, due to 
cancellations between positive and negative weights (see the last part on section ^). 

5.3 Comparison with data 

In this section, we would like to confront UMEPS with experimental data. Event generator 
predictions were obtained with the settings of Tune A2 [29] . The results should of course not 
be regarded as final statement, since changes in the perturbative physics in event generators 
in principle request a full re-tuning. The intention of this section is to investigate if after 
including matrix-element information, hard-scale features are closer to measurements, and 
to assess the changes in underlying event description. All plots were produced with RIVET 
[30]. We apologise if the selection of experimental measurements seems biased. 

In Figure |{J we show jet rates in W-boson production at the ATLAS [31]. We find 
an improved description after including up to two additional jets, and little differences 
between UMEPS and CKKW-L. 
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Figure 8: Transverse momentum of the third hardest jet, for pure QCD dijct production in pp 
collisions at Eqm = 7000 GeV, when merging up to two additional partons. Jets were defined 
with the /c^-algorithm, with kx,min = 10 GeV. Multi-parton interactions and hadronisation were 
excluded. The lower insets show the deviation of merged results from default PYTHIA8, for three 
different p M s-values. Left panel: Results of the CKKW-L scheme. Right panel: Results of the 
UMEPS scheme. 
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Figure 9: Jet multiplicity in W-boson production, for three different merging scales, as measured 
by ATLAS [31]. Effects of multiple scatterings and hadronisation are included. 



The transverse momentum of the hardest jet is shown in Figure 10. Again, we find 
that the shape description is improved by the CKKW-L and UMEPS methods. CKKW-L 
shows merging scale variations at lower pj_-values, since a slightly different inclusive cross 
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Figure 10: Transverse momentum of the hardest jet in W-boson production, for three different 
merging scales, as measured by ATLAS [31]. Effects of multiple scatterings and hadronisation are 
included. 

section for a low merging scale leads to a slightly different normalisations. UMEPS on the 
other hand suffers from statistical fluctuations for a low merging scale value (p MS = 15 
GeV), while the curves for p M s = 30 GeV and p MS = 45 completely overlap. Note that the 
p± spectrum of UMEPS is a little softer than CKKW-L, in accordance with Figure ||. 
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Figure 11: Azimuthal distance between the two hardest jets in W-boson production, for three dif- 
ferent merging scales, as measured by ATLAS [31]. Effects of multiple scatterings and hadronisation 
are included. 

It is interesting to investigate when tree-level matrix element merging schemes produce 
large uncertainties. Figure [Tl] shows the azimuthal distance A</>i2 between the two hardest 
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jets. The parton shower alone cannot describe the peak at ir. If the merging scale is low, 
the two-jet matrix element will give the dominant contribution in the peak region. High 
merging scales will increase the influence of the shower, thus degrading the description 
of the peak. Thus, this observable carries major merging scale variations, and provides 
excellent guidance for future improvements. 
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Figure 12: Jet shapes in QCD events, for three different merging scales, in two p± bins, as 
measured by ATLAS [32]. Effects of multiple scatterings and hadronisation are included. 
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Figure 13: Sum of transverse momenta of charged particles in QCD events, for three different 
merging scales, in the transverse and toward region, as measured by ATLAS [33]. Effects of multiple 
scatterings and hadronisation are included. 

Before moving to the discussion section, we would like to investigate how jet shapes 
at the LHC are changed by the inclusion of additional matrix elements in the pure QCD 
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case. The ATLAS analysis [32] found that the differential jet shape for relatively low-p_|_ 
jets with pj_ < 160 GeV depends crucially on the modelling of the underlying event. In 
Figure [l2|, we show the default PYTHIA8 and UMEPS results for two p± bins. Although 
far from perfect, the difference between the pure shower and merged results are similar to 
what a minor change in a s would give 6 . We are confident that the prescription for adding 
MPI (see section [4.1| ) does indeed mean that the underlying-event modelling of PYTHIA8 is 
only marginally perturbed by the inclusion of additional jets. This is supported by Figure 



13, which shows the sum of charged-particle transverse momenta in region close to the 



leading track (i.e. the toward region) and perpendicular to the leading track (transverse 
region) [33] . These are typical minimum-bias observables especially designed to investigate 
the underlying event, and PYTHIA8 tune A2 includes this data in the tuning procedure. It 
is reassuring that the inclusion of two additional jets through the UMEPS scheme did not 
invalidate this tuning. 



6. Discussion 



Before concluding this letter, we would like to make some comments to put this method 
into perspective. 

Relation to LoopSim 

Even though not completely obvious at first sight, UMEPS was heavily influenced by the 
LoopSim method [34]. This method as introduced to tame order-by-order large logarithmic 
enhancements by combining matrix elements with different jet multiplicities in a unitary 
way. The combination is done by joining all combinations to integrate over 1, . . . , n partons 
in the ME event S+ n , and also allowing an integration over hard process particles. In 
figure [I], the final state gluons are candidates for integration (or, in the terminology of 
LoopSim, looping), and the W— boson can be looped as well. With such a procedure, 
enhancements due to collinear W— boson radiation off a dijet state can be compensated 7 . 

Apart from major technical differences, one interesting difference is that in the LoopSim 
method, higher fixed-order corrections are approximated by multiple loopings, whereas in 
UMEPS, an all-order expression is included in the S +n+ i state before one parton is looped. 
Integrating multiple emissions is only necessary for /? MS — unordered sequences of splittings, 
which are not considered in LoopSim. Furthermore, in UMEPS, only integrations of QCD 
splittings are performed, while LoopSim includes loopings of W— boson radiation. We 
postpone the inclusion W— boson clusterings in UMEPS until a full electroweak shower is 
available in PYTHIA8. It would clearly be interesting to combine the Sudakov resummation 
in UMEPS with the multiple loopings done in LoopSim. To arrive at a better description 
of S + q configurations, one could e.g. take al contributions from the looped S + ± state, al? 
contributions from the twice-looped S+2, and all higher orders from Sudakov-reweighted 

6 This of course does not mean that we will to use this data for tuning, but that a s -choices for different 
tunes can have a comparable effect on p(r). 

7 At very large transverse momenta, the effect of multiple soft /collinear electroweak bosons becomes 
important. An appraisal of high-px observables in WZ- production has recently been reported [35]. 
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thrice-integrated S+3 events. The way the inclusive cross section is maintained in such a 
procedure will be less obvious than in the case of LoopSim or UMEPS. 



Merging scale dependence 

In the original CKKW-L algorithm it is evident that the dependence on the merging scale 
is absent to the accuracy of the PS. This means that for any observable, leading logarithmic 
terms on the form a™L 2n , where L = In pp/ p MS , are correctly cancelled to all orders. For 
a shower 8 which in addition is correct to next-to-leading logarithmic accuracy, also terms 
on the form a™L 2n_1 are correctly cancelled. 

What we have accomplished with the UMEPS method is that the total inclusive Born- 
level cross section is almost completely independent of the merging scale. In addition, 



if we look at the master formula in eq. ( 4.15 ), it is clear that for any phase space point 



(f) n with n resolved partons, the inclusive cross section, integrating all contributions from 
higher parton multiplicities, is simply given by B n . Although B n includes no-emission 
probabilities calculated by the shower, it has no dependence on the merging scale, and 
hence, all inclusive n-jet cross sections are independent of the merging scale. Since exclusive 
n-jet cross sections are the difference between the n-jet and n+l-jet inclusive cross sections, 
also these are independent of the merging scale. 

The independence is, of course, not exact for any real observable. A jet algorithm 
will not cluster an n + 1-jet state back to the precise n-parton phase space point as 
would the mapping of the parton shower, or symbolically for a general observable O, 
+n+ij)B n +i 7^ 0(S +n j) J s B n+ i^ n . However, as long as the observable is collinear- 
and infrared-safe this difference will not have any logarithmic enhancements, and as long 
as the n-jet state is well above the merging scale we can take this scale to be arbitrarily 
small, without changing the n-jet cross section. 

Also, as we have noted before, there are some n-parton states, such as the one in 
figure [l], which do not have an underlying n — 1-parton state reachable with a reconstructed 
PS emission. Unless the PS is amended with VF-strahlung splittings, such contributions 
will always give a small dependence. However, we have found these to be numerically very 
small in the cases we have investigated. 

Events with negative weight 

Contrary to the standard CKKW(-L) algorithms our new UMEPS procedure will produce 
negatively weighted events. There has in the past been a great reluctance in the experi- 
mental community towards using generators with negative event weights. Mostly this has 
been a question about problems in handling the statistics and that it seems wasteful to 
spend a huge number of CPU cycles to do a full detector simulation on an event, which in 
the end will be cancelled by another event with a negative weight. However, the acceptance 
for negative weights have increased, and today most experiments are using programs such 
as MC@NLO [36-39], which do produce a fair amount of events with negative weight. 



s Note that the shower in PYTHIA8 which is used in simulations in this paper has not been formally 
proven to be NLL-correct. 
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Clearly UMEPS is more wasteful than CKKW-L in this respect, and the number 
of events that need to be analysed to get the same statistics is more than doubled. In 
fact it can be shown that the variance in the event weights, when calculating the no- 
emission probability for the zero-jet case in CKKW-L with the Sudakov-veto algorithm, is 
proportional to IIs +0 — n| . The variance for UMEPS, where the corresponding factor is 
calculated by reclustering one-jet states multiplied by a no-emission probability, is of the 
form 



Thus, for small merging scales (giving small no-emission probabilities), UMEPS becomes 
very inefficient as compared to CKKW-L. 

We believe that the benefits of UMEPS outweigh this drawback. Also we note that the 
algorithm works in a way such that all events will either have zero weight or a weight of 
order (±) unity. This is because the no-emission probabilities are generated by the Sudakov- 
veto algorithm and are therefore either zero or unity, while the PDF- and a s -reweighting 
typically is of 0(1). Had the no-emission probabilities been calculated analytically, they 
would be very small for small merging scales, and e.g. each single 0-jet ME event would 
have to be cancelled by large number of small-weight reclustered 1-jet events. This would 
be very inefficient if a CPU-heavy detector simulation would have to be run on each event. 

On the other hand, the Sudakov-veto algorithm causes some problem in the case full 
detector simulation is not used. In this case the computational bottle neck is typically in 
the ME generation of high jet multiplicities, and the problem is that most of these events 
are given zero weight and will be thrown away by the Sudakov-veto algorithm, especially 
for small merging scales. This can in principle be handled by a modification of the veto 
algorithm [40] where all events are kept but are given a small weight. 

7. Conclusions and Outlook 

In this article, we have presented a new method for tree-level matrix element merging called 
UMEPS. This method is heavily indebted to the CKKW-L prescription, but explicitly keeps 
the total inclusive cross section fixed. Since it builds on the implementation of CKKW-L 
in PYTHIA8, all developments of CKKW-L are immediately available to UMEPS. This for 
example includes improvements for BSM processes [41] and multiple pre-defined merging 
scales. 

The UMEPS scheme uses an add-subtract prescription inspired by parton shower uni- 
tarity to combine the improved description of observable shapes of CKKW-L with a fixed 
total inclusive cross section. This means that significantly lower the merging scale values 
are possible, which allows for controlled improvement of low-scale features of the parton 
shower. Tuning efforts will be subject of a future article. When confronted with data, 
UMEPS and CKKW-L perform equally well. 

UMEPS is an ideal candidate for further improvements, since the lowest-multiplicity 
cross section is not reweighted, making replacements with the full NLO or NNLO results 




(6.1) 
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possible. We have successfully implemented an NLO extension, and will present it in a 
separate publication [42]. 

Finally, while finishing this article, it came to our attention that a very similar approach 
has been developed in parallel by Platzer [43], which further describes the extension of an 
inclusive-cross-section preserving merging method to NLO accuracy. 
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